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Abstract. We prove the homogcnization of the Dirichlet problem for fully nonlinear elliptic operators 
with periodic oscillation in the operator and of the boundary condition for a general class of smooth 
bounded domains. This extends the previous results of Barles and Mironescu [I] in half spaces. We 
show that homogcnization holds despite a possible lack of continuity in the homogenized boundary 
data. The proof is based on a comparison principle with partial Dirichlet boundary data which is of 
independent interest. 



1. Introduction 

We consider the homogenization of the following Dirichlet problem set in a bounded domain in 51 C K" 
with smooth boundary, 

F E (D 2 u e , x,-) = for xeQ 

(1.1) 

u £ (x) = g(x, f ) for x £ 90, 

with Holder continuous boundary data g £ C a (IR™ x K n ) for some a £ (0, 1). The operator (M, x, y) — s- 
F S {M, x, y) is assumed to be uniformly elliptic, and both the operator and the boundary condition y) 
are assumed to be Z n -pcriodic in the y variable. See Section 12.31 for the precise assumptions we make 
on the operators. 

The homogenization of the problem (jl.ip when g(x,y) = g(x) is independent of the fast variable is 
somewhat classical at this point and was done by Evans |121ll3j , the ideas originated in the unpublished 
paper of Lions, Papanicolaou and Varadhan [3D] on the periodic homogenization of Hamilton- Jacobi 
equations. There are too many works investigating the homogcnization of the interior operators in various 
settings for us to list all of them here, we simply mention that more recently the interior homogenization 
has been shown in the very general setting of stationary ergodic random media by Caffarclli, Souganidis 
and Wang in jS]. 

As far as oscillating boundary data is concerned there is less literature. For divergence form operators 
the case of oscillating co-normal Neumann data is in the classical book of Bensoussan, Papanicolaou 
and Lions [5] . For the case of non co- normal Neumann data and fully nonlinear operators there are a 
number of recent works. The papers of Arisawa [T] and Barles, Da Lio, Lions and Souganidis [5] identify 
the cell problem and show the homogenization in half-space type domains. The general domain case for 
the Neumann problem was solved by Choi and Kim in [9] (see also [10]). They require that the operator 
F e homogenizes to a rotationally invariant operator, in this situation the homogenized boundary data is 
continuous which guarantees the uniqueness of the limiting problem and thus stability of the e-problem 
as e — > 0. This issue is explained in more detail in [9]. Our work is based on the ideas in 9, 10j. Indeed, 
one of the novelties of this article is to show the uniqueness and stability result for a general class of 
operators even though the homogenized boundary data may be discontinuous [17j . In particular we do 
not need to put any special assumptions on the operators besides the uniform cllipticity. We believe 
that similar ideas should apply to the case of oscillatory Neumann data. 

For the oscillating Dirichlet data case there are even less references. For linear divergence form 
equations the homogenization in general domains was shown recently by Gerard- Varet and Masmoudi 
in [H1IT5] where they are able to consider linear systems, see also [T5]. The work of Barles and Mironescu 
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[4] deals with fully nonlinear non-divergence form elliptic operators when Q is a half-space. Barles and 
Mironescu identify the cell problem and the attendant difficulties in solving it under quite general 
assumptions on the operators. We will be concerned with extending this result to general domains. In 
order to clarify the main issues that arise from the general domain rather than from the cell problem 
we will avoid giving the most general assumptions on the operators. We expect that our results should 
also hold in the generality at which [4] show the cell problem can be solved, but this will be addressed 
in future work. 

At least on the surface a difficulty of this problem, especially in contrast to the corresponding os- 
cillating Neumann data problem, is a lack of equicontinuity for the u e . The optimal estimate for the 
continuity of the u e up to the boundary cannot be any better than the continuity of g(x, — ) which has 
oscillations of size 1 at arbitrarily small scales in the limit e — > 0. In particular the homogenization of 
(jl.ip . if it occurs, must be happening outside of a 'boundary layer' of width o(l). Keeping this in mind, 
let us call F to be the homogenized operator on the interior given by the results of [8j[T2l[l3]. Our goal 
is to show that the solutions u £ of the e-problem converge locally uniformly in to u solving an 
equation of the form, 



As we will see, the homogenized boundary condition g will depend not only on the function g, but also 
on v x , the inner normal to fi at x and on the F s . 

Let us describe heuristically the approach taken in [4] to identifying the homogenized boundary 
condition. The analysis of the problem proceeds by blowing up about points x £ dtt. As in [4], 
the introduction of the localizations, 



at points xq £ d£l in the limit as e is taken to leads us formally to the following "cell problem" set in 
the half space P(x ,v) := {y ■ (y - x ) ■ v > 0} where v £ S^ 1 \x\ = 1}, 



Some appropriate assumptions on the F s and g will yield that F is uniformly elliptic and that both 
F and ij) are Z"-periodic in y. Then by a Liouville type result the limit v{y + Rv) as R — > oo will 
be a constant ~fl(yp,F,v), we will identify this as the homogenized boundary condition. Then, at least 
formally, one expects that \u e {x + eRv) — v(Rv)\ —> as e — > 0. In particular, if we pull away from 
the boundary sufficiently far near xo, we can recover the interior averaging. This suggests that our 
homogenized boundary condition in (|1.2[) should be g(x) = Jl(g{x, •), F x , v x ) where v x is the inner unit 
normal to f2 at x and F x is an appropriate blow up limit of F e . 

The main difficulty that one has to deal with in extending to general domains is the rational boundary 
points x £ dfl such that the interior normal v x to f2 at x is in KZ™. Exactly when v £ RZ™ the cell 
problem (|1.3|) docs not yield a unique constant ~p. This corresponds to an a-priori lack of uniform control 
in v of the rate of convergence of v(y + Rv) — > ~jl{y). Some stability of this rate of convergence in v 
is necessary in order to show that the behavior of the cell problem governs the homogenization near 
irrational boundary points. 

It turns out that v(x + Rv) — > Jl(v) uniformly in v £ 5 n_1 is equivalent to Ji{ip, F, ■) having a unique 
continuous extension from the irrational directions S n ~ 1 \ IRZ™ to all of S 1 ™ -1 . When F is a rotationally 
invariant operator one can show using the methods of [9, 10] that such a continuous extension exists. For 
general F the homogenized boundary condition ~p is expected to have no such continuous extension [T7] , 
so we are left to consider weaker continuity conditions which still suffice to prove the homogenization. 




(1.2) 



v e {y) 



u £ {x + ey), 




(1.3) 
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One of the main contributions of this work is showing that the weak stability which one has 'intrinsi- 
cally' is sufficient for the homogenization to hold. We find this stability by carefully analyzing the rates 
of convergence for the cell problem in terms of the discrepancy function, sec Sections 12.21 and [51 

Let us call T the set of rational boundary points, 

r(O) := {x €dVt:i/ x e RZ n }. 

In light of the above discussion, the best continuity result we expect for general F e is, 

~g : dft — > K is continuous at x £ dil \ T, (1-4) 

and it turns out that d£l \ T is also exactly where one is able to exploit the connection between the cell 
problem and the e-problcm to show that any subsequential locally uniform limit of the u e satisfies the 
homogenized boundary data. In order to show that all subsequential limits are the same the fundamental 
question becomes, when is there uniqueness for the problem: 

F(D 2 u, x)=0 for x £ Q 

(1.5) 

u(x)=g(x) for x £ df2\r(0). 

We are not aware of any existing results in this direction, it is worth noting that this problem is related 
to the open question of whether there exists a boundary version of the Alexandroff-Bakclman-Pucci 
estimate. We provide a resolution of this question when r(fi) has sufficiently small Hausdorff dimension 
in the following Theorem, see Section [4] for a precise explication of all the assumptions, 

Theorem 1.1. [Theorem 14.2] Let F(M,y) be uniformly elliptic with constants (A, A) and satisfy ap- 
propriate continuity assumptions in (M,y). Suppose that f2 has the strict 7 exterior cone condition for 
some < 7 < 1. There exists /?o (n, A, A, 7) > max{0, 4(n — 1) — 1} such that if 

(3 > dim-u r(fi) > 0, 

then (jl.5p has a unique bounded solution u. 

Remark 1.2. The constraint on the Hausdorff dimension of r(f2) allows us to construct a sequence of 
super solution barriers <f>j which are large on r(f2), but converge to zero on compact subsets of Q. Then 
for any bounded solution v of (|1.5|) one can use the classical comparison principle to show v <u + cf>j, 
and letting j — > 00 will yield the result. 

Now we will give a statement of our main result. The results are not stated in their full generality in 
order that it be understandable at this point in the article. For the exact statements one can refer to 
the main text, in particular the main result is in Theorem 15. 31 

Theorem 1.3. (Solution of the Cell Problem and Homogenization) 

(i) [Lemmas 13.11 [3~31 and [3T4"] Let v XQM solve (jl.3|) in P(v,Xq), there exists a constant ~fl(ip , F,v) such 
that for all v £ S n ~ x \WL n , 

sup sup \v XoM (y + Rv) — Jt(tp, F, v)\ — > as R^-oo. 

io£i" y&P(x ,v) 

Moreover fx is continuous in all of its arguments at irrational directions v £ S*" -1 \ KZ™ (for the 
precise meaning of this refer to the referenced Lemmas), 
(ii) [Theorem 15.3] Let u £ be the solutions of the e-problem (|1.1[) and fl a bounded domain with C 2 
boundary satisfying the condition of Theorem \l.l[ Then the u 6 converge locally uniformly in f2 to 
the unique solution u of (|1.5j) with, 



sO) = v{g{x, ■),F x ,u x ). 
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Remark 1.4. In particular the homogenization result holds whenever the set of rational boundary points 
r is countable for general uniformly elliptic operators F e satisfying the assumptions of Section 12.31 
For example, by the inverse function theorem, this includes the classical assumptions of [H] for the 
homogenization of co- normal Neumann data, that dfl has no flat parts in the sense that D T v x (the 
gradient in the tangential directions at x) has rank n — 1. 

Remark 1.5. Here F x is a rescaling of the F e at x £ dft. Under the assumptions we will put on the F E , the 
operators F e themselves will not actually depend on e and one can simply take F X (M, y) = F e (M, x, y). 
For a precise characterization of the F x see Lemma 15.11 

Remark 1.6. It turns out that Jl(ijj, F, v ) is just the average of ip over a unit cell when F is linear and 
v is an irrational direction, see Section 13.11 for the proof. Since there is no dependence on the operator 
in the linear case one may wonder whether ~p is always just the average even when F is nonlinear. We 
show that in fact this is not the case in Section l3~Tl generically when F is nonlinear ~p{ip, F, v) depends 
on F and is not linear in the ip argument. 

Remark 1.7. We do not consider the case of divergence form operators here, for example, 



In this case the rescaled operator F x in the cell problem has the form, 

F x (D 2 u, Du, y) = -Tr(A(ai, y)D 2 u) - div y A(x, y) ■ Du. 

We believe that the only major difference to handle this kind of operator lies in the Localization Lemma 
12.71 As in [4] the requirement is that there exists a supersolution w x in the domain P(0, v x ) with 
w x (y + Rv x ) — > oo as R — > oo. In order for all of the continuity results to hold there will need to be a 
uniform lower bound on the growth of w x for x G dfl. For example this will hold for the above equation 
with the function w x (y) = y ■ v x as long as div y A(x, y) ■ v x > 0. 

1.1. Outline of the Paper. In Section [5] we discuss various notations and previous results which will 
be used throughout the paper. For readers well versed in the field of homogenization of uniformly elliptic 
equations most of the material in Section [2] except for Section 12.21 should be familiar. In particular we 
will discuss in detail the properties of Z" periodic functions restricted to hyperplanes with rational and 
irrational normal directions. These results are essentially refinements of the classical equidistribution 
theorem of Wcyl |5T] on irrational rotations of the torus. Then we describe precisely the assumptions 
on the differential operators we will consider and the associated comparison, regularity and interior 
homogenization results. 

In Section [3] we analyze the cell problem (JT73J) ■ We prove Theorem 11.31 part (i) and show properties 
of the homogenized boundary condition ~p in particular the continuity properties discussed above. We 
show that the convergence of the solutions of the cell problem v Xo>u (y + Rv) to JI is sufficiently uniform 
at irrational directions to show the claimed continuity. We also give an example which demonstrates the 
nonlinearity of the homogenized boundary condition when F is nonlinear. 

In Section [5] we use the results of the previous section to show the homogenization. In particular we 
use the stability of the uniform convergence v X0:U {y + Rv) — > Jl at irrational directions to show that the 
upper and lower half-relaxed limits of the u £ are respectively sub and supersolutions of (|1.5|) . Then we 
use the result of Theorem 1 1.1 1 to conclude that the upper and lower half- relaxed limits are the same and 
so the homogenization holds. 

1.2. Acknowledgements. I would like to thank my advisor Inwon Kim for reading this article and 
providing many helpful comments. I would also like to thank the anonymous referees for their suggestions 
and comments which helped to improve the clarity of the paper in various places. 
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2. Set Up 

2.1. Notations. Let us make a note of the notational conventions we will follow in this work. The 
numbers < A < A will always refer to the ellipticity constants. The meaning of the term ellipticity 
constants will be explained in Section [231 The number a G (0, 1) will almost always refer to the Holder 
continuity of whatever Dirichlet boundary data is under consideration. Let C M™ the associated 
Holder spaces are, 

C a (n) := (v : n -> K : U\\c°(n) ■= sup \^(x)\ + |V|c«(n) < 
where \ip\c a is the Holder semi- norm defined by, 

, ,, \^{x)-i>(y)\ 

\lp\c° := SUP ; . 

xjty&a F - V\ 

We also define the oscillation of a real valued function ib on K n over a set E C R™, 

osc ib := sup i/> — inf ib. 

e E e 

Constants denoted by C and c will denote universal constants, they depend only on A, A and the dimen- 
sion n, they may change value from line to line and sometimes even within the same line. Let A and B 
be two quantities, we write 

A < B to mean A < CB where C is a universal constant. 

If we want to emphasize the dependency of a constant on a quantity A we will write C(A). 

2.2. Rational and Irrational Directions. We define the following classes of normal directions with 
respect to the Z™ lattice, 

Definition 2.1. (i) v G S 71 ' 1 is called a rational direction if v G WE a . 
(ii) v G S n ~ x is called an irrational direction if it is not a rational direction. 

Given a Z™ periodic function ip we we will be interested in the properties of tp\ H , the restriction of 
ip to the hypcrplanc H = {x G W 1 : (x — xq) ■ v = 0}. When ^ is an irrational direction the distribution 
of the values taken by does not change too much under translation of H (changing xq). This 'fact' 
will be very important to understanding the homogenization of the boundary data in half spaces. 

We will first give a heuristic description of the issue at hand. Let us think of ip as a function on the 
unit periodicity cell [0, 1)™. Then we may consider the parts of the unit periodicity cell which are cut 
by the hyperplane H 

H per = {t£ [0, 1)" :T + zeH for some z G Z"} = |J (H + z) n [0, 1)". 

For hyperplanes H whose normal direction is rational the union defining H per will be equivalent to a 
finite union. As a result of this the restriction ip\ H ([0, 1)™) may be only a non-dense subset of all the 
values taken by %b. Moreover, this subset will be highly variable under changing xq- For example this 
issue is immediately evident in R 2 when we consider the periodic function, 

ip(x) = tp(xi) = sin27ra;i, 

and the hyperplanes H = {xi = a}. 

When the direction is irrational H per will be 'uniformly dense' (in a sense to be made precise) in 
[0, l) n . This is exactly equivalent to the well known irrational rotation of the torus example in the case 
7i = 2. In this case changing xq will not change the overall distribution of the values taken by ib\ H . 
To be slightly more precise, this is the situation in which homogenization of the oscillating boundary 
condition ip(x/e) will hold. 



6 



WILLIAM M. FELDMAN 



Now we proceed to give a precise description of the above heuristics in a way that will be useful to 
us in the future. Parts of the following discussion are borrowed from [lOj . The reason we repeat it here 
is we will need to make use of a more refined version of the results outlined there. 

Definition 2.2. A bounded sequence of real numbers (xj)j2. 1 is said to be equidistributed in an interval 
[a, b] if for any [c, d] C [a, b] we have, 

|{xi, ...,x n } H [c,d] | d-c 
hm = . 

n-*ao n b — a 

For x £ K let [x] be the largest integer less than or equal to x. 

Definition 2.3. A sequence is equidistributed mod 1 if the sequence (xj — is equidistributed in 

[0,1]. 

Theorem 2.1. (Wcyl's equidistribution theorem |21| ) If x is irrational then (jx)j e n is equidistributed 
mod 1. 

In order to make quantitative estimates we introduce the notion of the discrepancy. The following 
definition is from the book [18] via [10]. For a subset E C [0,1], a natural number TV and a sequence 
(xj)jzN C [0, 1] call A(E; N) the number of elements of the sequence Xj for 1 < j < N contained in E. 

Definition 2.4. For any x <G [0, 1] let Xj = jx — [jx] and define the discrepancy function, 



D N {x) 



sup 

B=[o,6)C[0,l] 



A(E;N) 



N 



\E\ 



Then by Weyl's equidistribution theorem Dn(x) — > as N — > oo for each x irrational. Now let us 
replace Dn{x) by a function equivalent up to constants which is continuous at every irrational. We 
define the modified discrepancy function as in |18) . 



D* N = D* n {x-l,...,x n ) = sup 

0<a<l 



A([0,a);N) 



N 



Then we have the following properties for D* N {x), 

Lemma 2.2. (Properties of D* N ) 

(i) (Theorem 1.3 in \18^) The discrepancies Dn and D* N are equivalent up to constants, 



D* N <D N < 2D* N . 



(ii) (Theorem l.J± in tlSjj ) Let x\ < X2 < • ■ • < xn be in [0, 1). Then their discrepancy D* N is given by, 



D 



N 



1 

2N 



max 
i=l,...,N 



2i-l 



2N 



(Hi) Define the modified discrepancy function for all x G M 

D* N (x) := D* N {x - [x],2x - [2s], ...,Nx- [Nx]), 

then this function is continuous in a neighborhood of every irrational i£l\Q. 

Proof. The proofs of the first two parts can be found in the book [15] . We prove part (iii). Given an 
irrational x and an integer N > let us call 



5 = min{|ja; - kx\ : 1 < j < k < N} > 0. 
Let a : {1, N} — > {!,..., N} be the permutation which orders the sequence x — [x], . 



,Nx-[Nx], that 



is. 



a{l)x - [a{l)x] < a{2)x - [a(2)x] < ■ ■ ■ < a{N)x - [a(N)x] 
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Then for \y — x\ < ^ we have that the above ordering is preserved for y, in fact for each 1 < j < N, 
o(j)y < o{j)x + <j{j)\y -x\ < a(j)x + 5/2 

< a(j + l)x - 5/2 < <j(j + l)x - a{j + l)\y - x\ 

< <r{j + i)y. 

Then from part (ii) we have that for \y — x\ < -Aj the discrepancy is has the form, 



2N 

Dtr(y) = h max 

N\y) 2iV i=i,...,JV 



a(i)y- [a(i)y] - ^j^- 



which, as a maximum of finitely many continuous functions which are continuous in a neighborhood of 
x, is also continuous in a neighborhood of x. □ 

Now let v = (yi, v n ) G 5 n_1 be a direction, and let 1 < i < n be the largest index corresponding 
to the largest component of v, that is, 

i = max{l < j < n : \vj\ = \v\£a°}. 

Let H u = {x £ M™ : x ■ v = 0} be the hyperplane through the origin orthogonal to v. For 1 < j < n 
define rrij(v) to be the slope of H u in the plane span{ei, ej}, 

= -mj(u)\ui\ + = -mj{v)\v\^ + \uj\. (2.1) 

The first equality indicates that one of the vectors mj{v)ei ± ej is in the hyperplane H u . The second 
equality makes manifest that rrij(u) thus defined are continuous in v and rrij{y) £ [0, 1]. Now for every 
e > and some 7 € (0, 1) we define, 

(J U (N):=2 min Dt m (mAv)) for N > 1. (2.2) 

l<j<n L J 

For irrational directions v G S'"~ 1 \Z n at least one of the rrij(v) will be irrational. As a result, cj v (N) —> 
as TV — > 00 when v is an irrational direction. Furthermore, for ./V sufficiently large, the min in (|2.2[) 
above will actually be the same as the min over j such that rrij(v) G K \ Q. Then for each such N > 1 
fixed tu.(N) will be continuous in a S"' 1 " 1 neighborhood of v with size depending on TV (it is a minimum 
of finitely many compositions of continuous functions). 

Now we state a lemma which quantifies our heuristic argument that hypcrplancs with irrational 
normals are uniformly dense modulo Z™ in [0, 1)™. The lemma is from [TU], but we state it in a modified 
form which will be useful to us. Due to the modifications we will also present the proof. 

Lemma 2.3. (Lemma 2.7 in [TO]) For v G S n ~ x and X{) G K™ we define H(xq) = {x G K™ : (x — Xq) ■ v = 
0}. Let uj v : TLyi — > M+ be as defined in (|2.2p . then the following hold: 
(i) There exists a dimensional constant C = C(n) > such that the following is true: for any x G 
H(xq) and any N > 1 there exists y G K™ such that 

\x~y\<C(n)N with y-x eZ n , 

and 

dist(y,if(.T )) < U) V {N). 
(ii) For any 5 > 0, there exists z G Z™ such that, 

distt z,H(x )) < inf loJN) +5. 

N>0 

Proof. Let v be any direction in S"" -1 and N > 1. Without loss we may translate so that x = and 
assume that f„ = |z^=o > l/^/n. Let j such that lo v (N) = f ^(mjfv)). Recall that one of the vectors 
mj(y)e n ± ej is in the hyperplane H = -ff(O), let us assume it is mj(v)e n + ej, the other case will 
work out the same. Now from the Lemma l2.2l we know that the the discrepancy D^(mj(v)) < uj v (N), 
and therefore every interval [a,b) C [0,1) with length at least oj v (N) contains at least one of the 
km,j(v) — [kmj(y)] with k < N. 
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Let s = s„e„, with s n G [0,1), then by the above argument there exists 1 < k < N such that 
\km,j{v) — \kmj{v)] — s n \ < oj v (N), and therefore, 

\k(mj(v)e n + ej) - s| mod Z" < u) v (N). 

In particular there exists y <E H u such that \y — s\ mod Z" < uj v (N) and |?/| < \/2N. 

Now let us consider s € [0, 1)™ arbitrary. Then let us take 

, s ■ v . 

s = e„ so that s — s € H. 

v n 

Now, by the above arguments, there exists y 1 € H such that \y' — s'\ mod Z™ < uj v (N) and \y'\ < \/2N. 
Then let y = y' + (s — s') which is still in H and \y\ < \y'\ + \s — s'\ < and, moreover, 

\y - s\ mod Z™ = \y' - s'\ mod Z™ < uj v {N). 

Now let s = xo mod Z™ , then there exists y € H as above with 

|y| < 2VnA^ and |y - x | mod IP = \y - s\ mod Z" < ^(iV). 

Now let y = y + (s — y mod Z ra ), then we have |y| < (2y/n + 1)N and 

y mod Z" = s = x mod Z™ and dist(y, ff) < ^(iV). 

This completes the proof of part (i). 

The proof of part (ii) is very similar. Let 6 > and N sufficiently large so that lo v (N) < inf jy>o u>v(N) + 
5. Let s = —xq and by the arguments above there exists 1 < k < N such that 

\k(mj(v)e n + ej) + x \ mod Z" < U) V (N) < mi^uj, y (N) + 5. 

Since xo + k{m,j{v)e n + ej) is in H(xo) this completes the proof of (ii). □ 

2.3. Assumptions on the Operators. Now we give the technical assumptions on the differential 
operators under which we can solve the cell problem. The operator F which arises in the cell problem 
is obtained as a scaling limit of the general operators F e being homogenized over. This connection 
will be made more explicit later in Lemma 15.11 We will work in the class of fully nonlinear uniformly 
elliptic operators without gradient dependence. This is mostly for convenience and we believe that our 
framework can be extended to the gradient dependent case by using the arguments of fj] . 

Let Ai n be the class of n x n symmetric matrices with real entries. Then the differential operator 
F : A4 n x R" — > K will be assumed to satisfy the following properties: 

(Fl) (Lipschitz continuity) F is locally Lipschitz continuous in Af" x R", and moreover there exists 
C > so that for all x, y £ W and M, N € M n , 

\F(M,z)-F(N,y)\ < C{\z - y\{l + \\M\\ + \\N\\) + \\M - N\\). 

(F2) (Uniform Ellipticity) There exists A > A > such that for all jel™, 

XTt(N) < F(M, y) - F(M + N,y) < ATr(iV) for all M, N e M n such that N > 0. 

(F3) (Periodicity) For every zePje K" and M e M n we have, 

F(M,y + z) = F(M,y). 

(F4) (Homogeneity) For alH > 0, M € M n and y £ R™ we have F(tM, y) = tF(M, y) and, in particular 
F(0,y) = 0. ' 

A wide class of operators of this form come from the theory of optimal control of diffusion processes, 
for example the Hamilton- Jacobi-Bellman operators, 

F(M, y) = - sup Tv(A a {y)M ) where A < A a (y) < A V a e A. 
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In much greater generality there are Isaacs operators coming from differential games which have the 
form, 

F(M, y) = - inf sup Tr(A afs (y)M) where A < A af! {y) <AV(a,/))eix6. 

In fact, all operators satisfying (F1)-(F4) can be written as an Isaacs operator, the proof is a basic 
exercise. 

The operator F £ is assumed to satisfy (Fl) with the local Lipschitz continuity in both the x and y 
variables, the uniform cllipticity (F2) and (F4) the F £ arc Z" periodic in y. The most general form of 
such operators is, 

F e (M,x,y) = f(x,y) — inf s\ipTr(A a ^(x,y)M) where A < A afs < A V (a, /3) G A x B. (2.3) 

Here the A a, ° have uniformly bounded Lipschitz norm, are periodic in their y variable, and satisfy the 
uniform cllipticity condition 

A < A af) < A V {a, /3) e A x B. 

The function / is in C 0,1 (R™ x R") and is Z™ periodic in its second variable. 

The standard notion of weak solution for the equation F(D 2 u, y) = in a domain 57 C R™ is the 
viscosity solution. The usual reference on the theory of viscosity solutions is [IT] . 

Definition 2.5. (Viscosity Solutions) 

(i) u : 57 —> R upper semi-continuous is a subsolution of F(D 2 u, y) < if for any <fi G C 2 (57) such that 
u — 4> has a local max at yo € Q, 

F(D 2 4>(y ),y Q )<0. 

(ii) u : ft — > R lower semi-continuous is a supersolution of F(D 2 u,y) > if for any (f> £ C 2 (fl) such 
that u — (f) has a local min at y G f2, 

(iii) u : 51 — >• R continuous is a solution of F(D 2 u, y) = if it is both a subsolution and a supersolution. 

The Dirichlct boundary data can be defined in the classical sense for viscosity solutions. We say that 
u is a subsolution of F(D 2 u, y) = in f2 with Dirichlet data i/j : <9f2 R if u is a subsolution and 

lim supu(z) < ip(y) for all y £ dfl. 

Supersolutions are defined analogously and solutions satisfy both the sub solution and the super solution 
property. 

One very useful property of viscosity solutions is their stability under uniform convergence. Let 
Uj : £1 — > R with j £ N for some domain (possible unbounded) 57 C R". Let us defined the upper and 
lower half-relaxed limit operations, 

lim sup* Uj{x) = lim sup{uj(y) : y £ Q and \y — x\ + j 1 < s} 

e— >0 

and 

lim inf* uAx) = lim inf{u,(y) : y G f2 and |y — x| + j^ 1 < e}. 

e— >0 

One can also have Uj : Oj — > R with the Qj — > 51 in the Hausdorff topology on compact subsets of R™ . 
One can easily check that if the upper and lower half-relaxed limits agree then the v? converge locally 
uniformly to their common value. In any case we have the following stability result for viscosity sub and 
supersolutions under the half-relaxed limit operations, see for more details, 
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Lemma 2.4. (Stability of viscosity solutions) Let Fj(M,y) : M. n x R" — > R be a sequence of operators 
satisfying (Fl) and (F2) which converge to an operator F(M, y) locally uniformly in M and uniformly in 
y G R™. Let CI be a domain in R" and let Uj : CI — ^ R be a bounded above (below) sequence of subsolutions 
(super solutions) of, 

F ] {D 2 u J ,y) <0 in CI ( Fj(D 2 Uj,y) > in CI ). 
Define u* the upper half-relaxed limit of the u 3 (resp. w* the lower half-relaxed limit), then we have, 

F(D 2 u*,y)<0 m CI ( F(D 2 u*, y) > in CI). 
Then we state the interior homogenization result in terms of the half-relaxed limits: 

Theorem 2.5. (Interior homogenization f8\ \ 12\ \T3 \/ ) Let F e (M,x,y) satisfying (Fl) with the local Lip s- 
chitz continuity in both the x and y variables, (F2) and F £ are Z" periodic in y. Let CI a bounded domain 
in R™ and u e satisfying, 

F e {D 2 u e ,x, < in CI 

(2.4) 

u e {x) < M on dCl 

for some M > 0. Then there exists an operator F(M,x) such that u* = lim sup* u e is a subsolution of, 

F(D 2 u*,x)<0 in CI 



u*(x)<M on dCl. 

The analogous result holds for supersolutions as well. 



(2.5) 



2.4. Maximal Operators and Comparison. Now wc recall the Pucci maximal operators associated 
with the class of uniformly elliptic operators with constants A, A. The basic results about the Pucci 
maximal operators can be found in the book [7]. For M G M. n we can always decompose M = M + — M_ 
with M± > and M + M_ = 0. We then define, 

V + (M) = ATriU+ - ATrA/_ and V~(M) = ATrM + - ATrM_. (2.6) 

Then from (F2) we may derive for M, N <E M n and x € R™, 

-V + (M -N) < F(M, y) — F(N, y) < -V~{M-N). (2.7) 

Note that by the inequality (|2.7|) any viscosity solution of F(D 2 u, x) = will satisfy, in the viscosity 
sense, 

-V+(D 2 u) < and - T~(D 2 u) > 0. 

Moreover, (|2.7|) implies that the difference w — u — v of a classical supersolution u and a classical 
subsolution v is itself a supersolution of the equation, 

-T-(D 2 w) > 0. 

This is also true of only semicontinuous viscosity sub and supersolutions using the method of sup and 
inf-convolutions originally used for this purpose by Jensen [16| . We state the Theorem in a slightly 
different form closer to what appears in [7], 

Theorem 2.6. (Comparison) Let CI be a domain in R' 1 not necessarily bounded. Let F satisfy (Fl) and 
(F2). Suppose that u and v satisfy, in the viscosity sense, 

F(D 2 u,y) < F(D 2 v,y) in CI. 

Then w = u — v is a subsolution of, 

-T + (D 2 w) < in CI. 
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In bounded domains the comparison principle for the Dirichlet problem follows from Theorem 12.61 in 
a standard way. A comparison principle for solutions of uniformly elliptic equations in half spaces with 
sublincar growth at infinity then follows from the following localization lemma which will be useful to 
us later, 

Lemma 2.7. (Localization) Let 5, e > and L > 1 and suppose that v is a viscosity solution of 

' -V+(D 2 v) < in P(0,e„)nQ L 

< v(y)<5 on dP(0,e n )nQ L (2.8) 

v(y) < L 1_£ on dQ L r\P, 

where Ql is the cylindrical region, 

Ql = {W\ < L} x {0 < x n < L} with x' = x — x n e n . 

Then we have, 

v(y) <S + 2±nL- e in Q 1 . 

Remark 2.8. Since the Pucci maximal operators are translation and rotation invariant the same result 
holds for the analogous problem in the half space P(xq, v) for any v € S n ~ x and xq € R". 

Proof. Consider the following barrier cj>, 

4>{x) = L-^Ux'l 2 - 2A( n - l)x 2 n ) + 6 + (2±(n - 1) + l)L~ e x n . 

It is straightforward to check that tfi(x) > v(x) on dQn and that is a smooth supersolution of 
— r P + (D 2 (f>) > 2A(n — l)L _1_e in x n > 0. From the definition of viscosity solution v < <f> in Ql- 
Therefore since 

4>{x) <6 + 2{±(n-l) + l)L- £ in Q u 

we get the result. 

□ 

A straightforward application of the localization lemma is the following comparison principle for 
bounded solutions of the half-space problem (|1.3[) : 

Lemma 2.9. Let ipi,ip2 £ C a (M"). Let v\ be a bounded upper semi- continuous subsolution of (|1.3|) 
with Dirichlet data ipx, and V2 be a bounded lower semi- continuous supersolution of (|1.3[) with Dirichlet 
data "02 ■ Then, 

sup (v 1 (y) - v 2 (y))+ < sup (0i (y) - ip2(y))+- 

yeP yedP 

2.5. Regularity Results. The Pucci maximal operators govern the worst possible behavior for solutions 
of equations satisfying the ellipticity condition (F2). In addition to the role they play in the comparison 
result, Theorem 12 .61 they are useful because regularity results hold uniformly in the ellipticity class 
(A, A). In particular we have the following result, for more details see the book [7]: 

Lemma 2.10. (Interior Holder estimates) Assume (F2) and (F4)- Let xo <E R n and v be a continuous 
viscosity solution of 

-V+{D 2 v)<0 and - V~(D 2 v) > 0. 
in B r (xo). For every a € (0, 1) there exists C = C(A, A, n,a) > such that, 

sup 'y-f <Cr-° sup v(x). 

x,yeB r/2 (x ) F D\ a xeB r (x ) 

By using barriers one can prove the following Holder estimate up to the boundary, 



12 



WILLIAM M. FELDMAN 



Lemma 2.11. (Boundary Holder estimates) Let g G C a (M. n ) for some a G (0, 1) and suppose that u is 
a viscosity solution of 

F(D 2 u) = in Bf = {\x\ < 1, x n > 0} 

(2.9) 

u = g{x) on T = dB+ n {x n = 0}. 

Then we have the following estimate up to the boundary, 

\\u\\ < C(n, A, A)(|| ff || C a (r) + sup M). 

° y- D i/2> B + 

3. The Cell Problem 

Now we begin discussing the actual homogenization problem. First we will analyze the cell problem 
and identify the homogenized boundary condition. We will work in the following half-space type domains 
in R N , 

P(v, x ) ={i£l":(i - xq) ■ v > 0}. 

Let ip : R n -> R be C' a (R n ) for some a G (0,1) and periodic with respect to the Z™ lattice. For 
convenience in the future we will sometimes refer to this class of functions as Cp er (M. n ). From here on 
a will always refer to the exponent of Holder continuity of the Dirichlet boundary data for the problem 
under consideration. Now we consider solving the following Dirichlet problem in P(u, Xq), which we will 
call the cell problem, 

F(D 2 v XOtV ,y)=0 in P{y,x ) 

(3-1) 

v Xo , v {y) = ip(y) on dP(i/,x ). 
A standard application of Perron's method in combination with the Localization Lemma 12.71 and the 
Comparison Principlc l2.9l give a unique bounded solution to the cell problem. Standard barrier arguments 
imply that the solution given by Perron's method achieves the boundary data continuously. 

Now let us define the upper and lower homogenized boundary conditions arising from the cell problems 
(EH), 

fJ,*(v, F, ip) := supjlim supv Xo ^(y + Ru) : y G dP(u, xq), xo G K™} 
(x4v,F,i)):=M{]imm£v XOiV (y + Rv):y€dP(v,xo), xq G M™}. v ' ; 

Then we call ~pL — + M*) to be the average of n* and /U*. Notice that fj,*, /i* and ~p are invariant 
under translations of ip. In fact we could have equivalcntly defined them as envelopes over ip(- + r) 
for t G [0, 1)™ rather than over the location of the hyperplane dP(v,xo). From now on we will work 
mostly with the quantity ~p(?Jj,F,i>) but the upper and lower homogenized boundary conditions will be 
mentioned occasionally as well since we will see that they are relevant at the rational directions. 

Now we will show that ~p — fj,* — /i» when v is an irrational direction and that the convergence of 
Vx ,v{y + Rv) to ~p is uniform in xq and y as R — > oo. Furthermore we provide a precise estimate on 
the difference fi* — /i» when v is a rational direction in terms of the discrepancy in the form of u> u . This 
lemma is analogous to Theorem 2.2 in [4]. 

Lemma 3.1. (Identification of the homogenized boundary condition) Suppose that v Xa = v Xo ^ is the 
unique bounded solution of p.lj) set in P(l>,xq) where v is any direction in S n ~ x and xq G K n . Let 
Jl(ijj,F,v) defined as above then we have for all N > 1, 

sup sup \v Xo (y + Rv) — TTj < C(n, A, A)(|^|c» + osciP)((N/R) a +Lj„(N) a ) for R>CN. (3.3) 

x £R n yedP 

In particular, when v is an irrational direction v G S n ~ x \ WL n we have that //* = /i* = Jl and for all 
directions v G S n ~ 1 , 

fi* (v, F, tp) — /i* (u, F, if)) < C(n,A,A)(|V|c- +osc^) mf w„(JV)°. (3.4) 
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Remark 3.2. The fact that the limit of v Xo ,u{y + Rv) as R — > oo does not depend on the location of the 
hyperplane dP (i.e. on xo) for irrational directions is essential for homogenization to hold. As mentioned 
before this is hcuristically because restrictions of ip to irrational hypcrplanes see all the values of tp. The 
Holder regularity of ip is not essential, but in general continuity seems to be necessary. 

Proof. We may assume without loss of generality due to invariance of the equation under adding a 
constant and scalar multiplication that \tp\c < 1 an d that Halloo = ^osc-0 < 1. 

1. First we show that i>o,i/(y + Rv) converges along subsequences locally uniformly in R" to constants. 
For now let us call v = vo, v , P = P(y, 0). By the comparison principle we have the uniform bounds, 

-1 < v < 1. 

Let i?fc -> oo be a subsequence such that v{R k v) converges as Rk — > oo to some constant /2. 

2. Due to Lemma \2. 31 for any p = p' + Rv G dP + Rv there exists z£ Z" such that, 

\z-p'\<CN and h := dist{z, dP) < U) V (N), 

Note that we have z ■ v = ±/i, we assume the + sign, the proof will be very similar in the other case. 
We define the translation of v, 

v(y) = v(y - z), 

which is a viscosity solution of, 

F(D 2 v,y) = F(D 2 d,y - z) = in P + hv 
v(y) = ijj(y — z) — ip(y) on dP + hv. 

Now for y £ dP + hv C P we have, 

Hy) - i>(y)\ < \v(y) - v(y - hv)\ + \^(y - hv) - ^(y)\ < Ch a . 

Therefore by comparison principle, 

sup \v(y) - v(y - z)\ < Ch a , 

y&P+hu 

and then in particular, by the interior C a estimates, 

\v(p + Rv) - v{Rv)\ < \v{p + Rv) -v(z + Rv)\ + \v(z + Rv) - v(Rv)\ 

< osc v{y) + Cuj v {N) a < C(uj v {N) a + (N/R) a ), 

\y-p\<CN 

where C is of course independent of the point p. Then using comparison, 

osc v(y + Rv) < osc v(y + Rv) < C{uj v {N) a + (N/R) a ). 

y£P y£dP 

Now for R k > R, since R k v e P + Rv, 

sup \v(y + Rv) - v(y + R k v)\ < C(^(N) a + (N/R) a ). 

yEP 

Letting R k — > oo in the above inequality we derive, 

sup \v , v (y + Rv) - fl\ <C{uj„{N) a + (N/R) a ). (3.5) 
yeP(0,v) 

3. Now we show that the estimate (|3.5p is independent of the spatial location of dP. Let xo G K n 
and v XQ be the corresponding solution of (|3.1[) in P(v,xq), and v as before be the solution of (|3.ip in 
P(v, 0). Let N > 1, due to Lemma [2~3l there exists z £ Z n such that 

dist((9P(>, x ) + z, dP(v, 0)) = h < u v {N). 

Then v(y) = v(y — z) solves the following problem in P(v, z), 

F{D 2 v,y) = F(D 2 v,y - z) = in P{v,z) 
v{y) = i/j(y - z) = ip(y) on dP(v, z). 
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Note that v trivially satisfies the estimate ()3.5[) . Let us assume that z <G P(u,xo), the other case will 
follow by a similar argument. For each y <G dP(v,z), we have that y — hv € dP(v,xo), so by the C a 
regularity up to the boundary for v Xo we get, 

\v xo (y) - V(y)| < \v xo (y) - v xo (y - hv)\ + |^(y - Zu/) - ^(y)| < Ch a . 

Therefore by the comparison principle Lemma 12. 9[ 

sup \v{y)-v XQ {y)\<Cuu{N) a , 

Combining the above estimate with estimate (|3.5[) for v, we get 

sup \v Xo>v {y + Ru)-ft\ <C(u v (N) a + {N/R) a ), (3.6) 

where the constant C depends only on the uniform ellipticity of F by way of the interior C Q estimates, 
and in particular not on the arbitrary point xq or N. Taking the suprcmum over Xq € K n and we get 
the desired estimate except with ft. From here we simply note that ft € [/j,*,/!*] so estimate (J3T5J) with 
ft replacing JI still implies (|3.4j) which in turn implies (|3.3j) with ~p (possibly at the cost of increasing the 
universal constant C(n, A, A)). 

□ 

We now list some properties of the homogenized boundary condition which can be derived in a 
straightforward way from the identification Lemma 13.11 

Lemma 3.3. (Properties of ~p) The following properties all hold for fi* and /i» as well unless otherwise 
specified. Recall that the distinction between the envelopes is only relevant at rational directions, 
(i) (Continuity with respect to F) Fix ip S Cp er (W n ), v € S n ~ 1 and let Fj(M,y) be a sequence of 

operators satisfying (Fl)-(F^) which converge to an operator F 00 (M,y) locally uniformly in M and 

uniformly in y, then, 

lim sup \p{il),F h v) -p(^,Foo,i/)| < C(n,A,A)(|V|o- +osct/>) inf u v (N) a . 

In particular this implies the continuity of fi with respect to locally uniform convergence of its F 
argument at all irrational directions v € S"^ 1 \ MZ". 
(ii) (Monotonicity) For any ipi,ip2 G Cp er (W l ) and F±,F2 satisfying (Fl)-(FJ i ) we have, 

tpi < ip 2 and F 2 < Fx implies ~p(ipi , Fi , v) < Ji(ip2, F*, v). 

For the following we fix v £ and an operator F satisfying (F1)-(F4), and we call JL(-,F,v) = /!(•) 

for simplicity of notation. 

(Hi) (Continuity with respect to ip) For ipi,ip2 € C2 er (W l ), 

|7Z(V>i)-mW>i)| < HVi-^IU- 

(iv) (Homogeneity) For any t > and ip € Cp er (W l ), 

Jiitip) = tJJ{ip). 
(v) (Constants) For ip G C^ er (R") and c e R, 

Ji(ip + c) = ~p(ip) + c. 
(vi) (Translation invariance) For any ip £ Cp er (W n ) and any t € [0, 1)™, 

-p{lp{-+T))=-p{lp). 

(vii) (Sub/ Super- additivity) If the operator F is convex (concave) then, for any ipi,ip2 € Cp er (R n ), 
fJ>*(ipi + 1P2) > + M*C02) (or Li*(ipi + 1P2) < V*(ipi) + H*(ip2))- 

In particular, if F is linear then ~p is linear as well at irrational directions. 
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Proof. All except (i) are easy applications of Lemma 13. II and the comparison principle. Let us therefore 
consider the situation in (i). As before we assume that \ip\c < 1 and that HV'lloo = ^oscip. Let 
~pj = ~pj(ip,Fj,Li) and let Vj for j G NU {°°} respectively solve the cell problems (|3.1|) for Fj in the 
domain P = P(0,v). By the method of half-relaxed limits [5] the Vj converge uniformly on compact 
subsets of P U dP to Voo. Let 5 > and let N sufficiently large so that ixj v (N) a < c6 + inf at>o ui v (N) a 
and then taking R > S- 1 / a C'N, we get from (pO|) . 

sup \vARv) --p.\<8+ Cuj„(N) a < 25 + C inf uj u {N) a . 
ieNu{oo} N >° 

Then fixing such an R, since Vj{Rv) — > v oc (Riy) as j — > oo, for j > J(<5, R) sufficiently large, 

l7*i - Pool < l«i W - Pi I + - «oc(^)| + IwooW - Pool < 5<5 + C inf ^(A) Q . 

Since <5 was arbitrary we conclude. 

□ 

Now let us fix our boundary data if) G C a (R") which is Z™ periodic and F(M,y) our operator 
satisfying (F1)-(F4). We will call ~fl(v) = ~p,{^),F,v) which is defined in Lemma [3.11 for all directions 
v G S' 1 ^ 1 \ KL n . We now consider the continuity properties of ~p with respect to the normal direction v. 
The goal is to show that ~p is continuous at every irrational direction. 

Lemma 3.4. Let v G S"^ 1 \ WL n . Then ~p satisfies the following continuity estimate for any N > 1 
and \v' — v\ < rj(i>,N) sufficiently small, 

|7T(i/) - Tl{v)\ < C(n, A, A)(Mo- + osc V)(iV^s \ v -j/\tfs+ u u (N) a ). 

Remark 3.5. In particular this estimate implies, since uj v (N) — > as N — > oo for irrational directions, 
that ~ft{y) is continuous at every irrational direction. At every rational direction we have an upper 
bound on the limiting oscillation of ~p in small neighborhoods of v. In fact using the methods of [51110] 
when the operator F(M,y) homogenizes to a rotationally invariant operator one can show that ~p has a 
unique continuous extension from S'™ _1 \ IRZ' 1 to all of S n . In general though, numerical experiments 
suggest that there are non rotationally invariant operators for which Ji will not extend continuously to 
the rational directions |17j . In that case the above continuity estimate may be the best possible. 

Proof. As before, without loss of generality (due to Lemma [33]) we may assume that IV'lc < 1 and that 
Halloo = iosc^. 

Fix v G S' 1 - 1 \RZ n and let v 1 G S^ 1 be any direction with \i/-v\ < r? < 1/2, let L, R with L» R 
and r\L << 1. The choices of r?, L and R will be made more explicit shortly. The idea is to explicitly 
compare the solutions of (|3.1j) in P(0, v) and P(0, v') and use the identification of the homogenized 
boundary condition Lemma 1331 

Let us call x' = x — (x ■ v)v, then we have, 

{x ■ v > \v' - v\L] C P(0, v') n {\x'\ < L}. 

Moreover, on {x ■ v = \v' — v\L] letting a(x') = —7^7- so that x' + a(x')v G dP(0,i/') we have that 
< \L\u' — v\ since r\ < 1/2 and by using the C a estimates up to the boundary, 

\v v {x) - v v i(x)\ < \v v (x) - v v {x')\ + \ip(x') - ip(x' + a{x')v)\ + \v u >(x' + a(x')u) - v v >{x)\ 
<C(n,\,A)L a \is' -v\ a . 

Then we apply the localized comparison estimate of Lemma 12.71 in the domain (P(0, v) + \v' — v\Lv) D 
P(0,i/) to get, 

\v v > -v v \< C(L a \v' - v\ a + £ ) in {\x'\ <R}x{0<x-is<R}. 
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Then by Lemma I3TI for every 1 < N < c(n)R, 

\m - W)\ < \v v (Rv) - m\ + C{L a W - v\ a + * ) + \v v ,{Rv) - -p{v')\ 

< C(L a \v' - u\ a + f + (f )« + w„(iV) Q + w„,(iV) Q ). (3.7) 
Now in order to make the first three terms on the right hand side above of comparable size we choose, 

R = N\v-v'\-^rs ail d L = \ v - v ' \~ (IXlfv+l) . 

So for any N > 1 choose 77 sufficiently small so that R > C(n)N, that is, 

< T) < C{n)- {2+a \ 

Then based on the continuity properties of uj. (N) with fixed N from Section 12.21 choose r\ smaller if 
necessary so that 

u V '{N) < 2u y (N) when \v-v'\ < r)(v,N). 
Then we get, by plugging in our choice of R, L for \v — v'\ < i](v, N) into (|3.7p . 

< C(n,\,A)(N^\is-i/\*k +u> v (N) a ), 

which was the desired result. □ 

3.1. The Nonlinearity of the Homogenized Boundary Condition. In the case when the operator 
F(M, y) = Tr(j4(?/)M) is linear the homogenized boundary condition ~p(ip) can be identified explicitly 
as the average of ip, 



7Zty>,*»= / iP{y)dy. (3.8) 
•/[o,i)« 

Note that in particular this value is independent of v £ \ WL n and of the operator F. The proof is 
a straightforward consequence of the Riesz Representation Theorem in combination with the properties 
given in Lemma 13.31 We state this in slightly greater generality in the following Lemma, 

Lemma 3.6. If F satisfying (Fl)-(FJ^) is concave then for all irrational directions v and all ip £ 

Cper 0^" ) • 



n(ip,F,v) > / if>(y)dy. 
J[o,i)« 

The opposite inequality holds in case F is convex. More precisely we argue that ifjt is linear in tp then 
it is just the average of ip over the unit periodicity cell. In particular this is the case when F is a linear 
operator. 

Proof. We will just show (|3.8[) . Then the result of the Lemma can be derived by noting that F concave 
(convex) with ^(0, x) = can be written as a supremum (infimum) of linear homogenous operators and 
using the monotonicity of Ji with respect to F. By Lemma [3~3l we see that JL : C a (T n ) — > M is linear and 
satisfies, 

\PW)\<W\U 

Therefore JI extends to a bounded linear operator (with norm 1) on the space C(T"), and so, by the 
Riesz Representation Theorem, there is a Radon measure a on T™ so that, 

m(V0 = / ip(y)do-{y). 



Then, from the translation invariance along with the fact that <r(l) = 1, we derive that a must be the 
uniform measure. □ 

Clearly this identification of fx, (|3.8|) . relies strongly on the linearity of the operator F. Nonetheless 
one may imagine, since JI is independent of the interior operator in the linear case, that this may carry 
over to the nonlinear case as well. This is not the case. We give a simple example which shows the effect 
of the nonlinearity: 
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Example 3.7. In dimension n = 2 let us define, for A > 1 and v, r\ orthogonal irrational directions, the 
following fully nonlinear and concave operator, 

F[D 2 u) = min{— Au, — u m — h.u vv \. 

We will take our boundary data to be ip(x) = cos(27rxi) and consider the solution u of the following cell 
problem, 

F{D 2 u) = in P{v) 

(3.9) 

u{x) = cos(27ra;i) on dP(v). 
Note that J, Q 1 % n ip{x)dx = would be the homogenized boundary condition associated with tjj if F were 
linear. In this case we can construct explicit solutions for both of the operators in the definition of F 
using separation of variables. In particular we consider, 

vi(x) = e~ x ' v cos(27T77i(a; • 77)) and V2(x) = e~ x ' v I A ' cos(27r?7i(x • 77)), 

which solve Q3.9[) with the interior operator replaced by Au and + Ku vv respectively. In particular 
V\ and V2 are both subsolutions of (|3.9p . Therefore, by comparison, we see that 

u{x) > max.{vi{x), ^(x)}. (3.10) 

Now let us consider the values of v\ and v-i for example on the hyperplane {x • v = A 1 / 4 }, here we have 

max{vi(x) 7 V2{x)} = e~ A ' [cos(27t?7i(.t • 77))]+ + e _A ' [cos(27r?7i(a; • fy))]-. 

Notice that this is the restriction to the hyperplane {x ■ v = A 1//4 } of the Z™ periodic function on IR™, 

i>'{x) = e- A_1/4 [cos(2Tr(xi - A 1/4 vi))}+ + e- Al/4 [cos(2tt(x 1 - A 1/4 Vl ))}.. 
Or in other words, 

ip'(x) = max{wi(x), V2(x)} < u(x) on dP(v, A 1 / 4 ^). 
Therefore using Lemma 13.61 identification of Ji for linear operators and the monotonicity of fl, 

/l(cos(27r-), F, v) > JL(ip', -A, v) = [ iP'(x)dx > i (e^ 17 * - e" Al/J ) > 0. 

J[o,i)" ^ v y 

As a result /I generically will not be just the average of the boundary data when F is nonlinear. Also 
note that in particular, since JI must be the average if it is linear in ip, it cannot be linear in ip. 

4. Comparison with Partial Boundary Data 

Recalling that we do not expect the homogenization to hold for at points on boundary with 

rational normal direction we are led to consider whether the comparison principle holds when the ordering 
on the boundary only holds on a 'large' subset. In particular one is led to consider the following situation, 
suppose that we have v and u bounded and respectively upper and lower semicontinuous in the closure 
of a bounded domain and a subset T C d£l such that, 

F(D 2 v, x) < F(D 2 u,x) in n 

(4-1) 

lim supw(y) < lim inf u(y) on <9f2 \ T. 
Then can we deduce that v < u in fl, or more generally is there an estimate of the form, 

sup (v - u) + < K |r|, 

xeK cesi 

where | • | is some measure of the size of T. Such an estimate, a form of the Alexandroff-Bakelman-Pucci 
estimate on the boundary, appears to be unknown even for linear F at least when | • | = "H™ -1 (-) the 
n— 1 dimensional Hausdorff measure. We are able to prove such an estimate only for | • | = H /3 °(-) with 
some /?o < n — 1, possibly quite small depending on the ratio A/ A. 
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In order to prove an estimate of the above form we will use the singular solution that was constructed 
in [2]. The following Theorem was proven there for even more general operators, 

Theorem 4.1. (Theorem 1 from [2]) For any 1 > 7 > let K 1 = {x : x n > (—1 + 7)|a;|}, there exists 



a 



unique constant /3o(n, A, A, 7) > max{j(fi — 1) — 1, 0} such that the problem, 



A > 

-V + (D 2 $)=0 in K, 

$ = on &KT 7 \{0} K ' 

has a positive solution $ G C(A' 7 \ {0}) which is — fto homogeneous, that is, 

$(x) = \x\~ Po ^(x/\x\). 

Note that by the Evans-Krylov Theorem $ is C 2 - a on the interior of for some a G (0, 1), see [7]. 
For rj G 5 n_1 let O v be any rotation sending 77 to e„ and then we call, 

$„(a:) = $(O n x). (4.3) 

In order to use the above barriers we will make the assumption that the domain f2 has the following 
exterior cone condition for some < 7 < 1, 

Definition 4.1. Let < 7 < 1, then a bounded domain il C K™ has the strict 7 exterior cone condition 
if there exists p > and 7' > 7 such that for every x G <9f2 there is a direction rj x with, 

{y:(y-x)- i lx < (-1 + j')\y - x\} n B p (x) C R" \ fi. 

Of course convex domains have the 7 exterior cone condition for every 7 < 1 since they have a supporting 
hyperplane at every boundary point. More generally any Lipschitz domain also has the strict 7 exterior 
cone condition for some 7 depending on the Lipschitz constant of dfl. 

We will also use the notion Hausdorff dimension below, so we will recall the definition. The d- 
dimcnsional Hausdorff measure is defined in the following way for a subset E C M. n , 

!oo 
VrJ : 3xj G l n s.t. £cll B{x J ,r j ) and r 3 < S 
;>i 

It is standard to check that W d (E) is decreasing in d and 7i d (E) = for d > n, this motivates us to 
define the Hausdorff dimension as, 

dim H E := inf {d > : U d {E) = 0}. 

Theorem 4.2. Let F satisfying (Fl) and (F2). Suppose that VL has the strict 7 exterior cone condition 
for some < 7 < 1 and let ^(n, A, A, 7) as above from Theorem 1 of J^. Let F C <9f2 such that 
- H' 3o (F) = 0, in particular this is the case if dim% F < /?o- Then if u and v bounded and respectively 
upper and lower semicontinuous on the closure of fl satisfy, 

F(D 2 v,x) < F(D 2 u,x) in n 

(4.4) 



lim supv(y) < lim inf u(y) on dfl \ T 



then u < v in fl. 



Remark 4.3. This Lemma is also true in the gradient dependent case as long as the operator is globally 
Lipschitz in the gradient argument. The only alteration is that the best estimate of /?o. In fact Theorem 
1 of [2] gives the existence of a singular solution in the gradient dependent case as well. 
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Proof. For simplicity we assume that the 7 exterior cone condition holds globally, that is p from the 
definition of the 7 exterior cone condition is equal to +00. A localization procedure using Lemma [2.71 
will cover the general case. Let w = u — v by the Comparison Theorem 12.61 it; is a bounded subsolution 
of 

-V + (D 2 w) < in n 

(4.5) 

lim sup w(y) < on dil \ V. 

Let $ be the singular solution from Theorem 1 of [2] corresponding to 7, then call m = min{$(0) : 8 6 
K 7 > n S" 1 ^ 1 } > and take M > m~ 1 \\w\\ 00 . Then from the definition of Hausdorff measure, T-L l3a (T) = 
implies that for each 8 > there exists B(xj,rj) for j e N, an open covering of T by balls, so that 
^i>o r j° ^ /M. For each Xj let rjj be from the 7 exterior cone condition for O and consider the 
following barrier, smooth in f2, 

M*)=5 + J2 Mtfi Vj (x -Xj). 

We have the following convergence of the (j>$ \ 0, 

0<(/>s{x)<28 for d(x,dn) > 8 1/p . 

Let us show that <f>s is a strict smooth supcrsolution on the interior of f2, since {y : (y — Xj) ■ r\j > 
(-1 + i)\y - x 3 \} C £l c for each j along with (|4~2j) , 

-V + {4>s) >J2- Mr j V+ ( D2 ^v } ( x ~ x o)) > 0- 
i>o 

Let us show that for each a; € dfl, 

lim sup (w — 4>s){y) < —8. 

This is true by the assumption for x € dtt \ T since <ps > 8 on <9f2. For x € T we have that x £ B(xj, rj) 
for some j > 0, so that, as above, for y G 17 with |y — ir| < d(x, dB(xj, rj)), 

Mv) >$ + MtjQru (y - x 3 ) >8 + Mr^mrJ >S+ HU 

This proves the ordering of the Dirichlet data. Then since cf>$ is smooth in the interior of f2 simply by 
the definition of viscosity solution we have that w < <ps m Letting 8—^0 gives the result. 

□ 

5. HOMOGENIZATION IN GENERAL DOMAINS 

Now we consider the homogenization of in general domains f2. We show that (jl.ll) homogenizes 
when Oft does not have too many rational boundary points in the sense of Theorem 14.21 Unlike in the 
situation when D, is a half space which was considered in [1] it is now necessary to make use of the 
continuity properties of the homogenized boundary condition ~p with respect to v G S" -1 . The first 
reason is in order to guarantee that the homogenized problem, 

F(D 2 u, x)=0 in CI 

(5.1) 

u(x) = g(x) on dfl. 



has comparison. Here F is given by Theorem 12.51 and we will identify ~g{x) by the results of the previous 
section. This is where we will use the result of the previous section on comparison with partial boundary 
data. 

In order to blow up at boundary points and use the results for the cell problem in a locally uniform 
way we will again need the stability of the rate of averaging at irrational directions. We make this precise 
in Lemma \5 . 2 1 which is based on Lemma 3.1 in [3], but first we make the following note on the behavior 
of F £ under the blow up rescaling: 



20 



WILLIAM M. FELDMAN 



Lemma 5.1. The following convergence holds: 

e 2 F e {e~ 2 M, x, y) -> F x ( M, y) as e -> 0, 

locally uniformly in M and uniformly in (x,y), with F x (M,y) satisfying assumptions (F1)-(F4). In 
particular we also have that if Xj — > x then FJ — > F x locally uniformly in M and uniformly in y. 

Proof. By our assumptions on F e in Section [^751 there exists a collection of matrices, XI < A a @(x, y) < 
AI, A a '@ have uniformly bounded Lipschitz norm and are periodic in their y variable, also there exists 
a function / <G C 0,1 (R™ x M. n ) periodic in it's second variable so that, 

F £ (M,x,y) = inf sup [-Tr(A^(x, y)M)] + f(x, y). 

This form of the operators makes the claim clear, and in this situation, 

F x {M,y) = inf sup [-Tr(A afi (x, y)M)] . (5.2) 

□ 

Lemma 5.2. (Behavior of u e away from the boundary layer) Suppose that Q is a bounded domain with 
<90 being C 2 . Let xq € <9f2 such that the interior normal v XQ is an irrational direction. For any 6 > 
there exists Rq and rj > such that for R > Rq and for any sequence x e — > x with \x — xq\ < r), 

lim sup \u E (x e + eRv Xo ) — g(x)\ < S, 

(in particular the quantity in the lim sup is defined for e sufficiently small) where the homogenized 
boundary data g is given by, 

g(x) :=Jt(g(x, ■),F x ,v x ). 

Proof. First we explain how to choose rj. Given <5 > 0, and v Xo an irrational direction, there exists 
N Q > 1 sufficiently large so that oj„ x (N ) < 5 x / a . Recalling that oj. is continuous at irrational directions 
and that v x is continuous in x since Q is assumed to be a C 2 domain, let r\ be sufficiently small so that 

\x - x \ < ?y implies u Va (N) < 6 1/a . 

Now fix x € dn be fixed as above with |cc — xo | < rj as above and for every e > define 

v £ (y) := u £ (x + ey). 

Let us call fi £ = e _1 (il — x). The v £ satisfy the following equation in the viscosity sense, 

!F e (e~ 2 D 2 v e ,x + ey,s~ 1 x + y) = for y e tt e 
(5.3) 
v E (y) = g(x + ey, e 1 x + y) for y £ dQ s . 

The f2 e — > P x = P{Q,v x ) in Hausdorff distance when restricted to compact sets. Now we use the 
compactness afforded by the periodic setting. Let us call 

[0,l) n 3t e ^e^x mod Z", 

i.e. e~ x x — t £ + z e with z e £ Z™. We may replace e~ 1 x in (|5.4j) above by t £ using the Z n periodicity of 
F £ and g. Then taking a subsequence of the e the r e —> r G [0, 1]™. Noting that — ||g||L°°(R"xK") < v e < 
||<?||l°°(r™xr™) w e may take the upper and lower relaxed limits of v s along the above subsequence, 

v*' t (jj) = lim sup*i> e (2:) and v*, T {y) = hm inf * i; e (2:). 

The upper and lower half-relaxed limits are defined for each y £ P x . Let us show that these are 
respectively sub and supcrsolutions of the limiting scale one equation, 

F x {D 2 v t ,ij + t) =0 for y E P x 

(5.4) 

Vr(y) = g{x,T + y) for y e dP x . 
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The interior sub/supersolution property follows from the stability of viscosity solutions under uniform 
convergence, Lemma liOl The condition on the boundary can be shown in the following way. Fix y £ dP x 
and S > 0. For e sufficiently small, from the local Hausdorff distance convergence of d£l £ — > dP x , there 
exists y' £ dfl e n B(y, S). Now for any z £ Q £ n B(y, S). by the estimates up to the boundary Lemma 

nn ' 

\v e {z) - g(x, r + y)\< C\\g(x + e-,r e + .)l|c-l* - y'\ a + \g(x + ey',T e + y) - g(x, r + y)\ 
<C(\\g\\ c «)(5 a +e a \y'\ a + \T £ -Tn. 
Then taking the supremum over such z and letting e — > along any subsequence where r 6 — > r, 
\v*' T (y)~ g{x,T + y)\ < lim sup \v £ (z) - g{x, r + y)\ < C\\g\\ c «S a . 

£ ^° ze^ e r\B(y,d) 

Letting 6 — > proves that v*' T achieves the boundary data in (|5.4[) . The same argument works for v*. T . 



Since the equation (|5.4[) has comparison we derive v*' T < tVr, of course the opposite inequality is 
always true so we get the local uniform convergence in P x along subsequences of v £ — > v T for some 
t £ [0, 1]™ depending on the subsequence. The point is that the limiting behavior of v T (y + Rv) will not 
depend on r £ [0, 1]™. 

Let us consider v T = v T (y — t), these satisfy, 

F x (D 2 v T ,y) = for y £ P x + (r • v x )v x 

(5.5) 

Vr{y) = g(x, y) for y £ dP x + (r • v x )v x . 
Then by Lemma [3TT1 for any N > 1, 

sup sup ^(y + iJO -pCffCa, 0,^,^)1 ^ ^(iV) a + (iV/i?) a . 
re[o,i]" ye-Px+(r-^)i/ x 

Given 6 > fixing N = No as in our choice of above and letting i?o > sufficiently large we get, 

sup sup \v T (y + Rv x ) - g(x)\ < S. 

t£[0,1]" y€dP x +(T-v x )v x 

Then for R > 2i?o, f] > sufficiently small by the continuity of v x so that \v x — v XQ \ < 1/2, and each 
r£ [0,1]", 

\u e {x + eRv Xo )-g{x)\ < \u £ (x + sRv Xo ) - v T (Rv Xo )\ + \v T {Rv XQ ) - g(x)\ < \u £ {x + eRv Xo )-v T (Rv)\+5, 
and infimizing over r and then taking lim sup on both sides, 

lim sup \u £ (x + £Rv XQ ) — ~g(x)\ < S + lim sup inf \u £ (x + eRv x „) — v T (Ru Xo )\ = S. 

e^O e->0 re[0,l] n 

Since the convergence of the u £ is locally uniform along subsequences we may replace x by any sequence 
x £ — > x above. □ 

From Lemma 15.21 we will be able to deduce that homogenization of the boundary condition holds 
near boundary points with irrational directions. In general one cannot blow up to get an estimate of 
the form in Lemma 15.21 near boundary points with rational normal directions. However we have the 
comparison principle with partial boundary data from Theorem 14.21 This Theorem gives a condition 
under which satisfying the Dirichlet data only at boundary points with irrational normal directions is 
enough to deduce uniqueness. 

Now let C M. n be a bounded domain with C 2 boundary which has a strict 7 exterior cone condition 
for some < 7 < 1. Recalling that v. : dfl — > S"" 1 is the inward unit normal vector field to dtt we call 
the set of rational boundary points, 

T = T(O) := {x £ dQ, : v x £ RZ n }. 
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Let us suppose that O has no flat boundary pieces in any rational direction of too large Hausdorff 
dimension, precisely we assume, 

■H Al (r) = o, 

where n — 1 > (3o(n, A, A, 7) > max{0, 4-(n — 1) — 1} is the singular solution exponent from Theorem 1 
in [2]. In particular this is always true of uniformly convex domains which have an exterior half space 
at every boundary point and for which T will be countable. More generally any bounded C 2 domain 
such that the set of rational boundary points where D T v x (the tangential gradient) is degenerate has 
Hausdorff dimension smaller than (3q will satisfy this condition. 

Theorem 5.3. Suppose that £1 satisfies the hypothesis given above. Letg(x) = ~fi(g(x,-),F x ,v x ), then 
~g : <9fi — > K is continuous at every point of dQ\T, and the solutions u e of (jl.ip converge locally uniformly 
in Q to the unique, in the sense of Theorem \4-2\ solution given by Perron's method u of 

F(D 2 u,x)=0 in n , . 

u(x) = g{x) on dtt \ F. ^ ' 

Proof. Let us first note the continuity of ~g at points of <9£1 \ T. Let Xj G 9f2 converging to x G dfl \ F, 
then ||ff(a;j, •) — g{x, < C\x n — x\ a , v Xn — > v x since dfl is C 2 and F x * (M, y) converges to F X (M, y) 

locally uniformly in M and uniformly in y as j — > 00. Then from the continuity of Jt noted in lemma [3751 
(i) and (iii) and in lemma HOI we have the convergence, 

V( x j) ='P{9{xj,-),F x ^v Xj ) -> -p(g(x, -),F x ,v x ) =g{x) as j -> 00. (5.7) 

This shows the claimed continuity of g. 

The main tool to show the homogenization is Lemma [5721 Fixing xq € dQ. \ S, v XQ the interior normal 
to dQ. at xo, 8 > and R > Ro(S) and 77 = rj(6) > from Lemma I5T21 we consider, 

u e (x) := u £ (x + sRvq) defined for x G — bRvq. 

Since the u e are uniformly bounded (by comparison with iH^H^) the upper and lower relaxed limits 
are defined and finite in Q, 

u* := lim sup* u e and u* := lim inf * u 6 . 

Note that for x £ the upper (and lower) relaxed limits of u 6 and u e (which we call u* and w* 
respectively) agree. The only dependence on R, the manifestation of the boundary layer for u e , is on 
dQ. By our choice of R and r\ from Lemma 15.21 and the continuity shown above in ()5.7|) there exists 
< if < 77 sufficiently small so that for x G dfl with \x — xo\ < if , 

u*{x) < g{xo) + 8 and u*(x) > g(xo) — 5, 

and by the interior homogenization result for x G we get the inequalities, 

F(D 2 u*,x) < and F(D 2 u*, x) > 0. 

Let 4>± be respectively super and subsolution barriers for the domain D, at xq. Then, as in the standard 
barrier argument, we consider 

<p+ + {g{xoj + 0), 



\x — Xq I , 

which satisfies (f>+ > u* on dQ. By comparison (f> + > u* in Q. In fact since <f>± can be chosen smooth 
this is even just by the definition of subsolution. In any case, since u* = u* in D, we get, 

lim sup u*{z) < lim sup u*(z) < g(xo) + S. 

Z—tXQ Z^Xq 

Making the same argument with 0_ and we get, 

lim inf u*(z) > g(xa) — 5. 

Z^Xq 
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Of course now the left hand side above is independent of 5, so we may take S — > to get, 

g(xo) < lim inf u*(z) < lim sup u*(z) < g(xo). 

This argument works for every xq € dQ \ T. Applying the comparison principle with partial boundary 
data Theorem 14. 2 1 for F to u* and u* we get for x € fi, 

u*(x) — u*(x), 

so that the u e converge locally uniformly in D, to u = u* = . 

□ 
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